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'■ Abstract: 

^ : 

We show that the Brodsky-Lepage evolution equation for the spin 3/2 baryon distribution 
Qh. amplitude is completely integrable and reduces to the three-particle XXXs=-i Heisenberg 
r-| ! spin chain. Trajectories of the anomalous dimensions are identified and calculated using the 
^1 expansion. Extending this result, we prove integrability of the evolution equations for 

' twist 3 quark-gluon operators in the large limit. 

^3: 
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A QCD description of hard exclusive processes invokes the concept of hadron distribu- 
tion amphtudes (DAs) given by matrix elements of nonlocal light-cone operators bewteen the 
vacuum and the hadron state. It was realized long ago that one-loop renormalization group 
(RG) equations for the leading twist meson distributions are diagonalized in the conformal 
basis: The mixing matrix for the corresponding local operators becomes diagonal and the 
anomalous dimensions coincide with those familiar from studies of the deep-inelastic scat- 
tering. The situation with baryon DAs is much more complicated. In this case one has to 
consider renormalization of three-quark operators of the type: 

B{a, 6, c) = eijkq\au)q^{hu)q^{cu). (1) 

Here g is a quark field and u is an auxiliary light-like vector = 0; taking the leading-twist 
part and insertion of gauge factors is implied. Conformal symmetry allows to eliminate 
all mixing with operators including total derivatives but is not sufficient to diagonalize the 
mixing matrix. In the usual expansion of baryon DAs in Appell polynomials only the lowest- 
order term (asymptotic DA) is renormalized multiplicatively while in general one is left with 
a nontrivial (A^ + 1) x (A^ + 1) mixing matrix, A^ being the number of derivatives, and has 
to diagonalize it explicitly order by order, see [@-^. The analytic structure of the spectrum 
was not known and further analytic results were not expected. 

Main result of this letter is that the 3-particle Schrodinger equation describing the renor- 
malization of spin 3/2 baryon operators is completely integrable, i.e. has an additional 
integral of motion. Our result is similar in spirit to the recent discovery [^,0 of integra- 
bility of the system of interacting reggeized gluons in QCD, but is obtained in a different 
physical context. In particular, the ERBL evolution equation for spin 3/2 baryon operators 
appears to be mathematically equivalent to the equation for the odderon trajectory, and the 
results obtained in the latter context can be adapted to unravel the spectrum of baryon 
operators. 

Extending this result, we prove integrability of the RG equations for twist 3 quark- gluon 
operators 

S^{a, b, c) = q{au) [iG^,^{bu) ± G^y{hu)'-f^]u''i q{cu), (2) 
T{a,b,c) = q{au)u^u^a^^Gyp{bu)V q{cu), (3) 

where F = {/,i75}, in the limit of large number of colors N^.. Such operators give rise to 
twist 3 nucleon parton distributions and have attracted considerable interest recently, see 
e.g. [§|10|- Further results will be presented elsewhere [p!T[| . 



To one-loop accuracy, the divergent part of the nonlocal operator $(01,02,03) where 
$ = 5, 5*^, T, has the form (l/e)if$(aj), and the explicit expression for the integral operator 
H is known for all cases under consideration P, pr2| , |l0| . An arbitrary local operator O with 
A^ covariant derivatives can be represented by the associated polynomial in three variables 
ipi^ai, a2, as) of degree A^ such that = ip{da, db, (9c)$(a, b, c)a,b,c^o where da = d/da etc. In 
order to find multiplicatively renormalizable local operators one has to solve the Schrodinger 
equation for the V'-functions, Hip = Sip, where H is easy to find if H is given. 

It proves convenient to define the integral transformation |T3| ip{ai) ipi^i) by 
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where k and si are the canonical dimension and spin projection of the i-th field, respectively: 
I = 3/2 , s = 1/2 for quarks (antiquarks) and / = 2 , s = 1 for gluons. We can reformu- 
late the above eigenvalue problem in terms of ip functions and it is easy to check that the 
corresponding Hamiltonian H coincides with the initial Hamiltonian for the nonlocal oper- 
ator, H = H. As a trivial consequence of topology of one-loop Feynman diagrams H has 
a two-particle structure, H = J2i,kHik- Conformal invariance implies that the two-particle 
Hamiltonians Hik commute with SL{2) generators J^'^ = I]i=i jf 



''^'^ where 



Jf = zA + {k + Si)/ 2, 



and are hermitean with respect to the scalar product 

(^l|V^2) =iJl{di,... ,(93)V^2(^1,... ,Z3] 



Zi=0 



(5) 



(6) 



Thus, the equation Hip = Sip decays into the set of eigenvalue problems on the subspaces of 
functions with fixed value of J3, J^ip = j^ip and annihilated by J_, J_iIj = 0. The advantage 
of the '^-representation' is that in this basis the latter condition is simply shift- invariance 
T3||. Therefore, eigenfunctions of two-particle Hamiltonians are given by simple powers 

(Zi 



Zk) instead of Jacobi polynomials in standard variables ||14 



The SL{2) invariance imposes stringent restrictions on the form of two-particle operators, 
so that only a few structures are allowed. One such structure corresponds to the 'vertex 
correction' involving the gluon field from (one of) the covariant derivatives (in Feynman 
gauge): 



a I 



ll){z'^^, Z2, Z^) -1p{z) 

ip{zi, Z21, Z3) - ipiz) 



(7) 



where z = {zi, Z2, z^,}, z"^ = Zia + z^a and a = 1 — a. Another structure originates from 
gluon exchange between quarks (or between a quark and a gluon): 



Hl^iPiz) 



2 jDa^y-—^iP{z^i,z^l,z,) 



(8) 



where Da = IlLi dai 6{1 — !]««)• 

Due to SL{2) invariance the two-particle hamiltonians must depend on the corresponding 



Casimir operators Lik = {Ji + Jk 



Jik{Jik + 1) only. This dependence can be easily 



reconstructed from the spectrum of Hit. Since the form of the eigenfunctions is known 
ipi = [zi — ZkY, it is straightforward to derive 



h:,^'"^ = 2 [ipij^k + 1) - ^(2)] , 

H^k^''^ = ^{-hk + 3/2) + ip{.hk + 1/2) - ip{?>) - ^(2), 



(9) 



where ipi^x) = r'(x)/r(a;). The superscripts (qq) and (qg) indicate quark-quark and quark- 
gluon operators, respectively. Similarly, we obtain 
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Hti"'^ = 2{J,u + 3/2)-i( J,, - l/2)-\ (10) 

We are now in a position to specify RG equations for the operators in (|l]-0) explicitly. 
One has to distinguish 3-quark operators belonging to (3/2, 0) and (1, 1/2) representations, 
which correspond to DAs for spin 3/2 and spin 1/2 baryons, respectively. We get 

H,/2 = Hli'^'^ + Hli"'^ + Hli'^'\ (11a) 
= ^3/2 - (l/2)i/rf^ - {l/2)Hlt'^. (lib) 

Omitting subleading in Nc terms, the quark-antiquark-gluon Hamiltonians are |[9|,p^ 

Hs^ = Hli''^ + Hlt'^^ - Hlt'\ (12a) 
Hs- = ifit"^) + Hli''^^ - Hli''^\ (12b) 
Ht = Hli'^'^ + Hli"''^ - Hlt'^ - Hli"''^ . (12c) 

The properly defined anomalous dimensions are given in terms of eigenvalues of the above 
operators including color factors and trivial contributions of self-energy insertions: 

73/2,1/2 (iV) = (1 + 1/iVc) ^3/2,1/2 (iV) + (3/2) Cf, 

75,T(iV) = £sAN) + (7/2) (13) 

where Cp = {N^ - l)/(2iVc). 

The operators Hs± are equivalent; hereafter we consider Hs+- The corrections to 

(|1^) and RG equations for 3-gluon operators involve additional structures [ill] ^^^^ 



be discussed here. 

We have been able to find integrals of motion Qi, [Hi,Qi] = 0, for all Hamiltonians 
in question with the exception of Hi/2- Explicit expressions for the conserved charges Qi 
present the main result of this letter: 

Q3/2 = z[Li2,Li3] = «9i(92<93Zi22;23^31, (14) 

9 1 

Qs+ = {Li2, L23} - -L23 - 2^12, (15) 
9 9 

Qt = {L12, L23} — 2^12 — 2^23, (16) 

where { , } stands for an anticommutator. Remarkably, -^3/2 coincides in our representation 
with the Hamiltonian of the XXXs=_i 3-particle Heisenberg spin chain. Hence the expres- 
sion in (p!4D for the conserved charge Q3/2 follows directly from the corresponding classical 
result, see also @]. 

To check that [Ht, Qt] = 0, [Hs, Qs] = we introduce a complete set of functions |13 



ip^k{zi, Zk] zi), n = 0, . . . , N where z, is a fixed pair of indices (say, (1,2)) and zi is the third 
variable, different from Zi,Zk. The functions ip^i^ are obtained by the 'hat '-transformation 
(I) of the set of polynomials of degree N: 
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N-n 
m=0 

Sk) N—n—m 

where CJ^ is the binomial coefficient and {a)k = T{a + k)/r{a). ipnizi, Zk) is defined such that 
^pn{zi, Zk) = {zi — Zk)"' is an eigenfunction of Hik and the factor Z„ is chosen by requiring 
that ip^k has unit norm. 

The functions ip"k are shift invariant and mutually orthogonal with respect to the scalar 
product (H). It is easy to check that Ja'i/'j'^ = (iV + 7/2)-^,'].. The operator Lik is diagonal 
in the basis 'ip^/^ while the other two Casimir operators are three-diagonal {{n\L\n') ^ for 
\n — Ti'l < 1 only). 

Consider matrix elements of the commutator 

A = [Qt, 4''^] = {Li2 - 9/4, [H^''\ L23]} 

where H^'^^ = £[^2'^^^ — H\2^^\ sandwitched between the -^"2 'states'. Since L23 is three- 
diagonal in this basis, the only nonzero elements are A^^n+i and An+i^n- Due to antiher- 
miticity of A it is sufficient to consider An^n+i- 

An,n+1 = {L23) {£n — £n+l){^n + — 9/2), 

(12) 

where Sn and are the eigenvalues of the operators Hj, and L12, respectively. Using 

explicit expressions for Hj. , L12 it is easy to derive that An,n+i = 2(L23)n,n+i('^^n — -^n+i), 
that is in operator form 

[QT,Hi}^^]=2[Li2,L2s]. (18) 

Similarly, we obtain that [Qt,h!^^'^] = 2 [1723,-^12] and, consequently, [Qt,Ht] = 0. The 
proof for operator is analogous. 

Once conserved charges are known, one can consider the eigenvalue problem for these 
charges instead of the Hamiltonians, which is simpler. For the Heisenberg spin chain, a 
detailed study exists due to Korchemsky The spectrum of Q3/2 is shown in Fig. la. 

For generic integer N there exist + 1 eigenvalues which come in pairs ±g. Note that for 
even Q3/2 has zero eigenvalue q = 0. The corresponding value of energy can be calculated 
exactly: 

^3/2(iV,g = 0) =4^(iV + 3)+47s-6, (19) 

see the dotted curve in Fig. lb. Eigenvalues of Q3/2 lie on trajectories which were found in 
using a 'semiclassical' expansion in the parameter h = N + 3: 



(0) 



q{N,k)/h' = Y.q^"'\k)/h"^, 

m 

1/^27, q^^\k) = -{k + l)/V3,... (20) 
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The q^"^\k) are polynomials of degree m; first eight of them are given in Eq. (5.14) in Ref. 
/c is a nonnegative integer which numerates the trajectory. The asymptotic expansion in ( 
is valid for g > only and the analytic continuation of the trajectory to g < can be obtained 
by using symmetry properties of the solutions 0] : 

q{N, k) ^ -q{N,N-k) 

Two examples of the trajectories with k = 2 and k = 7 are shown in Fig. la together with 
exact eigenvalues (crosses) calculated numerically. Note that the two asymptotic expansions 
— for positive and negative q — match reasonably well. Explicit analytic formulas for the 
trajectories in the g — >■ region are available from 

The low-lying eigenvalues of q{N,k) can be calculated to 0{1/N) accuracy from the 
equation |1TT| 

(q/N^) \nN- arg[T{l + iq/N^)] = -^(^ - 2A;) (21) 

which is valid for k — N/2 <^ InN. The lowest value of |g| for odd N is thus of order 

q/N^ = ±^(lniV + je)-' + 0(l/(lniV + 7^;)^). (22) 

The spectrum of -^3/2 is shown in Fig. lb. Exact eigenvalues obtained by explicit di- 
agonalization of the mixing matrix are shown by crosses. Since £^3/2(9) = ^3/2{—(l) all 
eigenvalues except for the ones for g = are double-degenerate. The energy eigenvalues lie 
on trajectories corresponding to the trajectories for q in Fig. la, and, similar to the latter, 
can be calculated using a 'semiclassical' expansion 

00 

r(0) _ 



E,/2{N,k)=e^'^ -Y.e^^\k)/h^, 

m=l 

=61n(A^ + 3) + 67£;-6-31n3,... (23) 



The polynomials e^™^^{k) are given in Eq. (6.5) of Ref. up to m = 7. The trajectories 
corresponding to = 2 and k = 7 are shown in Fig. lb by broken lines, whereas the solid 
curves correspond to the asymptotic expansion in (^) |jl6| . Note that the expansion diverges 
close to the 'deflection points' which occur at even integer N and with the energy given by 
Eq. (|19D . Convergence of the 1/h expansion is somewhat worse for the energy compared to 
the conserved charge q, but it can be improved systematically. Alternatively, one can derive 
analytic approximations for the conserved charge q{N, k) and £^3/2(0') applicable in the g ^ 
region, see 

Using (p2D one can derive an estimate for the lowest energy eigenvalue for odd A^: 

^3/2(iV) = 41niV - 6 + 47,, + j^^. (24) 

Numerically the difference to Eq.(^) is very small, compare exact eigenvalues with the 
dotted curve in Fig. lb, and is probably irrelevant for phenomenological applications. One 
has to bear in mind, however, that an approximation of taking into account operators with 
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the lowest anomalous dimension only for each is theoretically inconsistent since they 
belong to different trajectories. 



The anomalous dimensions of quark-gluon operators (plSf), (|T6]) can be studied along 



similar lines ||rT|. To the 0{\/N) accuracy the spectrum of low- lying energy eigenvalues is 



given in terms of eigenvalues of the corresponding integrals of motion as 

^(z/) = 2 In iV + 47£; - 5 + 27^e [V^(3/2 + iv)] (25) 

where 2z/| = qs^r — 3/2, and quantization conditions for the effective charges read, to the 
same accuracy 

TTTi 

z/rlnA^ + <l>i(z/) - $3(z/) = — , 

1 TTTl 

vs\nN + -{^i{v) + - ^^{v) = — (26) 

where n = 1,2,... and 

<l>i(z/) = (1/2) arg [J^^{'i/2 + iv, -3/2 + iz/, 1 + 2iv, 1)], 
<l>2(z/) = arg [^i{l/2 + iu, -1/2 + zz/, 1 + 2w, 1)], 

^^{v) = arg[T{'i/2 + iv)]. (27) 

These formulas are not applicable to the exact solutions with minimum anomalous dimen- 
sions, found in Refs. PJTO[], which correspond to imaginary v and have to be treated sepa- 
rately. We can show that these special solutions are separated from the rest of the spectrum 
by a finite gap. A detailed study is in progress. 

To summarize, we have shown that a few important three-particle evolution equations 
in QCD are exactly integrable, that is they possess nontrivial integrals of motion. This 
allows for a fairly complete description of the spectrum of anomalous dimension of baryon 
operators with spin 3/2, and similar techniques can be developed for other well. The 



eigenfunctions can also be studied |g,|TT[. We believe that the approach based on integrability 
may find many phenomenological applications to studies of higher-twist parton distributions 
in QCD. 
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FIG. 1. The spectrTim of eigenvalues for the conserved charge Qs/2 (a) and for the spin 
Hamiltonian -^3/2 (b), see text. 
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